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Abstract 

The XXZ model (f = 1 /2) in a transverse field on a double chain with a uniform long-range interaction among the z components of the spins 
is considered. The nearest-neighbour interactions are restricted to the components in the xy plane and to the spins within the same chain leg, 
such that the Hamiltonian is given by H = - ELl 'm L%i + £ A+l ) - h Em,«=l I^=i S z mJ S^ k - KfX=i L%i %n,j> where N 

is the number of sites of the lattice and m,n (m,n = 1,2) label the chain legs. The model is solved exactly by introducing the Jordan-Wigner 
and integral Gaussian transformations, which map the Hamiltonian in a non-interacting fermion system and corresponds to an extension of 
the model recently studied by the authors for a single chain. The equation of state is obtained in closed form, and the critical classical (at 
T > 0) and quantum (at T = 0) behaviours can be determined exactly. The quantum critical surface is determined in the space generated by the 
transverse field and interaction parameters, and the crossover lines separating the different critical regimes are also obtained. It is also shown 
that, differently from the results obtained for the single chain, the system can present multiple quantum transitions. 
© 2008 Elsevier B.V. All rights reserved. 
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The critical behaviour of the XXZ model with long-range 
interactions among the spins along the z direction, in a sin- 
gle chain, has been recently considered by the authors [ 1 ] . In 
that paper they present a comprehensive study of the quantum 
and classical transitions undergone by the model. Although the 
model presents first- and second- order quantum and classi- 
cal transitions, which is a very important feature, no multiple 
transitions are present. However, multiple quantum transitions 
can occur in magnetic systems as it has been shown in the ex- 
perimental results recently obtained for ZrZn2 [2], This result 
has been the main motivation for studying spin models which 
are exactly soluble and can present multiple quantum transi- 
tions. Therefore we will consider the XXZ model (s = 1/2) 
in a transverse field on a double chain with long-range uni- 
form interaction among z components of spins. The nearest- 
neighbour interactions are assumed isotropic, and restricted to 
the spin components in the XY plane and to the spins within 
the same chain leg. It should be noted that due to the isotropy 
of the nearest-neighbour interactions the total magnetization 
along the z direction is conserved. Assuming periodic boundary 



conditions the Hamiltonian is given by 
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where N is the number of sites of one leg and m ( m = 1, 2 ) 
labels the chain legs. By introducing the Jordan-Wigner trans- 
formation [3] and a Gaussian transformation, we can write the 
partition function as 
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where P = 1 /kgT , the effective Hamiltonian H is given by 
H = £ §J m (c] n j c mJ+l +4 J+1 c mJ )+£( $h + v^p7)4 jCmJ) 
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(3) 

and it has been used the fact that the long-range interaction 
commutes with H. Then, introducing the spatial Fourier trans- 
form, the Hamiltonian H , apart from a constant, can be written 
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in diagonal form with energy excitations E^ m = $J m cos(k m ) + 
P/z + y/lfilx. Then the partition function can be recast in the 
form 



1/2 
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n[l+exp(-£0] (4) 

m,k 



and, in the thermodynamic limit, can be evaluated by the saddle- 
point method. Following this procedure, the free energy func- 
tional can be written in terms of the magnetization, M z , and at 
T = 0, is given by 

f = h + 4IM z (M z + l)--(Ji sin 9] + / 2 sin(p 2 ) - - (cpi +(p 2 ) , 

(5) 

where M z is explicitly given by 
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with cp m =arccos(-(/z+4/M z )/7 m ). It should be noted that for 
Ji = 7 2 the results obtained so far reproduce the known ones for 
the single chain [1]. By minimizing the free energy functional, 
given in eqs (5) and (6), and the by using stability condition 
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we can determine the equation of state and obtain the phase 
diagram at T = 0. Whithout loss of generality, we will restrict 
the analysis of the phase diagram to the region J2/J1 ^ 1 , which 
presents the main features of the model. In order to charaterize 
the different quantum phases displayed by the model we present 
in Fig. 1 the magnetization as a function of the field for different 
values of the long-range interaction. As it can be seen, for / < 
0, the transitions are of second order since the magnetization 
is continuous at the transition points. On the hand, for / > 
the magnetization is discontinuous at the the transition points 
which is a signature of a first order transition. 

The quantum critical surface obtained is shown in Fig. 2 
where, as already mentioned, we have quantum phase transi- 
tions of first and second order. In the region / < 0, it can be 
shown that there are two critical surfaces of second order quan- 
tum phase transition. These critical surfaces are given explic- 
itly by the expressions h = J\— 2/arccos (—J1/J2) /n, which 
corresponds to the magnetization M z = l/27iarccos(— J1/J2) , 
and by h = / 2 — 21, which corresponds to the magnetization 
M z = 1 /2. These two surfaces intercept when J\ = / 2 as it can 
be seen in the diagram shown in Fig. 2. For / > the model 
presents three lines corresponding to the quantum phase transi- 
tions of first order, with a triple line which is also shown in Fig. 
2. In this case, it is not possible to obtain an explicit equation 
for the critical surface and it has been determined numerically. 

The interception of the critical surfaces of first and second 
orders occurs at / = 0, and corresponds to a bicritical line since 
it contains the known bicritical point present for J\ = J 2 [1,4], 
when the results reduce to the ones for the single chain. 

Moreover, it should be noted that, for a given set of param- 
eters, the system presents multiple transitions as a function of 
the field, as it can be seen in Fig. 2. In the second-order region, 



the multipe transitions can also be obtained from the explicit 
equation for the critical surface given above. 

Finally, we would like to point out that, for J\ — 7 2 , the 
quantum critical behaviour is identical to the one obtained for 
the model when we have a single chain [1,4], and as already 
mentioned we do not have multiple transitions. 
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Fig. 1 . Magnetization (at T = 0) as a function of the field for different values 
of the long-range interactions. 
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Fig. 2. Phase diagram of the model at T = and J2/J1 > 1. 



